F/B  20/11 


A RAY  THEORY  FOR  ELASTODYNAMIC  STRESS  INTENSITY  FACTORS. (U) 

MAR  7T  JO  ACHENBACH*  A K 0AUTESEN  N0001A-76-C-0063 


UNCLASSIFIED  NU-SML-TR-77-1  NL 


(/■  / ■ ■ - 

' — ^ J.  D./Achenbach  a«d  A.  K.^Gautesen 


The  Technological  Institute 
Northwestern  University  1 
Evanston,  Illinois  60201 


I 


ABSTRACT 

Elastodynamic  stress  intensity  factors  generated  by  the  interaction  of 
wave  motions  with  a crack  are  analyzed.  It  is  shown  that  in  an 
asymptotic  approximation,  which  is  valid  for  iTigh  frequencies,  the  stress 

| 

intensity  factors  at  the  edge  of  a crack  are  related  to  the  fields  of 
incident  rays  by  a matrix  of  stress  intensity  factor  coefficients,  which 
can  be  computed  from  canonical  solutions.  The  canonical  solutions  are 
provided  by  the  fields  describing  diffraction  by  a semi-infinite  crack 
of  plane  body  waves  and  plane  surface  waves,  which  are  incident  under  an 
arbitrary  angle  with  the  edge  of  the  crack.  Several  applications  of  the 
theory  are  presented.  For  cracks  of  finite  length,  the  contributions  due 
to  the  travelling  back  and  forth  of  rays  between  the  two  crack  tips  is 
taken  into  account  in  a simple  manner,  to  yield  results  which  are  in 
excellent  agreement  with  numerical  results  obtained  by  other  authors. 


1.  Introduction 


A crack  in  a solid  body  gives  rise  to  singular  stresses  at  the  edge  of 
the  crack.  In  a local  coordinate  system,  with  the  origin  on  the  edge  of 
the  crack,  a representative  stress  component  near  the  edge,  say  t may  be 
expressed  in  the  form 

Tp(r,p)  = (2rrr)'^  Ky  T?(0)  (1) 

Here  r and  g are  polar  coordinates  in  the  plane  normal  to  the  edge.  The 
geometrical  features  of  the  crack  geometry,  as  well  as  the  parameters 
describing  material  properties  and  loads,  enter  in  the  stress  intensity 
factor  Kj.  For  dynamic  problems  K^.  also  depends  on  the  time  t and  the 
circular  frequency  w. 

There  is  much  interest  in  stress  intensity  factors,  because  they  are 
relevant  to  fracture  mechanics  considerations.  Here  we  are  interested  in  elasto- 
dynamic  stress  intensity  factors  generated  by  the  interaction  of  wave  motions 
with  a crack.  We  present  some  general  results  which  are  valid  for  high  fre- 
quencies under  time-harmonic  excitations,  or  for  small  times  after  the  arrival 
of  wave  fronts  under  impact  loads. 

Ray  theory  provides  a very  useful  method  to  analyze  both  the  propagation 
of  high  frequency  waves  and  the  propagation  of  surfaces  of  discontinuity, 
see  e.g.  Refs. [11  - T3].  The  general  ideas  of  ray  theory  were  extended  by 
Keller  [41  to  analyze  the  far  field  in  diffraction  problems.  In  the  present 
paper  a ray  theory  for  stress  intensity  factors  is  established.  The  theory 
states  that  in  an  asymptotic  approximation. the  stress  intensity  factors  at 
a point  on  the  edge  of  a crack  are  related  to  the  field  of  the  incident  wave 
at  this  point  by  a matrix  of  stress  intensity  factor  coefficients,  which 
depends  only  on  the  direction  of  incidence,  the  frequency,  and  the  local 


physical  properties  of  the  solid  at  the  point  of  incidence.  The  incident 
wave  may  be  a body  wave  as  well  as  a surface  wave.  The  matrix  is  independent 
of  the  curvature  of  both  the  wavefront  of  the  incident  wave  and  the  edge  of 
the  crack,  provided  that  these  curvatures  are  of  smaller  orders  in  magnitude 
than  the  wavenumbers. 

In  Section  2,  an  edge  region  analysis  is  presented  for  incident  body 
waves.  Details  are  worked  out  for  an  incident  ray  of  longitudinal  motion. 

By  matching  a boundary  layer  solution  of  the  diffracted  field  in  the  vicinity 
of  the  edge  to  the  outer  solution  (the  geometrical  optics  field),  it  is  shown 
that  the  stress  intensity  factors  generated  by  the  incident  field  are  related 
in  a simple  manner  to  the  stress  intensity  factors  for  a plane  longitudinal 
wave  incident  upon  a semi-infinite  stress-free  crack.  The  modifications 
required  to  treat  incident  rays  of  transverse  motion  are  obvious,  and  they 
are  not  discussed  in  detail. 

For  incident  body  waves  the  construction  of  the  matrix  of  stress  intensity 
factor  coefficients  is  discussed  in  Section  3.  The  coefficients  are  listed 
in  Appendix  A.  In  Section  4 the  theory  is  extended  to  elastodynamic  stress 
intensity  factors  generated  by  surface  waves  which  are  incident  on  the  edge 
via  the  faces  of  the  crack. 

In  the  last  section  several  applications  of  the  theory  are  presented, 
and  the  results  are  compared  with  those  obtained  by  other  authors. 

2 . Edge  Region  Analysis 

In  this  section  a formal  proof  of  the  theory  for  incident  body  waves 
is  presented  by  a method  motivated  by  the  work  of  Buchal  and  Keller  r51. 

An  inner  solution  is  constructed  which  is  asymptotically  valid  in  a region 
near  the  edge.  The  leading  term  of  this  solution  is  then  matched  to  the 
leading  term  of  the  outer  solution,  which  is  the  geometrical  optics  part  of 


the  diffracted  field. 


Essential  to  the  development  of  this  section  is  the  "geometrical 
optics"  part  of  the  diffracted  field,  for  the  case  that  a plane  longitudinal 
wave  is  diffracted  by  a semi-infinite  crack,  whose  surfaces  are  free  of 
tractions.  The  geometry  is  shown  in  Fig.l.  Omitting,  here  and  in  the 
sequel,  the  factor  exp(-iUJt),  where  UJ  is  the  circular  frequency,  the 
incident  wave  is  represented  by 

vine  = £ exP<ikLP-£>  (2) 


In  Eq.(2),  we  have  k = ou/c  , where  c is  the  velocity  of  longitudinal 

Li  III  Li 

A 

waves:  cT  = r(\+2u)/o"l  , and  p is  a unit  vector  which  defines  both  the 

direction  of  propagation  and  of  displacement.  In  terms  of  the  angles  shown 
in  Fig.l  we  have 

(P1>P2»P3)  = (sin^L  coseL,  sin<JL  sin0L,  cosiJL)  (3) 

The  geometrical  optics  part  of  the  diffracted  field  has  been  worked 
out  in  a paper  by  Achenbach  and  Gautesen  [6].  The  results  may  be  written 
in  the  form 


* 


w 

~go 


= exp(ikLP3x3)  Wg0(Xl,x2) 


where 


(4) 


. 


V (xl’x2)  = Htnc  H(0-6L)  + ~L  H(e+6L"2n)  + ~T  H<6+V2n) 
where  H(  ) is  the  Heaviside  step  function.  The  relevant  angles  are 

defined  by 

CL  C0S^T  = CT  C0S<^L 

c^  sin0T  coseT  = cT  sin^  cos0L 

In  these  relations  c^,  = (p/o)2  is  the  velocity  of  transverse  waves. 
The  terms  in  Eq. (5)  are 


(5) 


(6) 

(7) 


mm  .«■«»»  »> 
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S inc  = P exp[ikL(Plx1+p2x2)] 

WL  = Rl  PL  exprikL(p^Xl+p2x2) ] 

T T 

WT  = Rt  d exp[lkT(p1x1+p2x2) ] 

where  and  R are  the  usual  reflection  coefficients  for  the  reflection 

of  a plane  longitudinal  wave  which  is  incident  on  a stress-free  surface 

under  an  angle  with  the  normal, and  krf  = u;/c  . We  have 

2 2 

sin20o  sin202-(cL/cT)  cos  202 

~ 2 2 
sin20Q  sin202+(cL/cT)  cos  202 

2(cl/ct)  cos202  sin20o 


rt  = 


2 2 

sin20o  sin0T+(cL/cT)  cos  20T 


(8) 

(9) 

(10) 


(ID 


(12) 


The  unit  vectors  appearing  in  Eqs.(8)  - (10)  are 
PL  = (PpP^P^)  = (Pi»‘P2’p3) 


(13) 


(sin0T  cos0T,  - sin0T  sineT,  cos 0 ) 
T Tv  , T, 


(14) 

(15) 


d = p x (i  x p ) / li  x p 

/V  ~ Z A/  *vZ  ~ 

Now  let  us  consider  an  incident  longitudinal  wave  whose  surfaces 
of  constant  phase  are  arbitrary  but  smooth.  The  wave  is  incident  on  a 
stress  free  crack  of  arbitrary  shape,  whose  edge  is  a smooth  curve.  The 
geometry  is  shown  in  Fig. 2.  We  consider  a point  0 on  the  edge,  and  we 
define  an  orthonormal  coordinate  system  (x^,x2>s).  Here  s is  arc  length 
along  the  edge,  with  positive  direction  such  that  the  propagation  vector 
of  the  incident  wave  makes  an  acute  angle,  0 (s),  with  the  tangent  to  the 

Li 

edge.  The  positive  axis  is  in  the  crack  surface,  and  in  the  direction 
of  the  principal  normal  to  the  edge.  The  x2~axis  is  in  the  direction  of 
the  binormal  to  the  edge,  and  such  that  the  coordinates  (x^,x2>s)  form  a 
right-hand  coordinate  system.  At  point  0 the  incident  wave  may  be  represented  by 
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u = A(x  x s)  yS  exp[ik  S(x  x ,s)l 
'»inc  1 z L l,  z J 

where  SCx^^^jS)  is  the  phase  function  which  defines  surfaces  of  constant 
phase  or  wavefronts.  The  rays  are  normal  to  the  wavefronts.  At  point  0 
the  direction  of  the  ray  is  defined  by  the  propagation  vector  p(s),  where 

A • 

P ( s ) = Vs(0,0,s)  = (sin0  cose  , sin0  sin©  , cos0  ) 

a/  L L L L L 

Thus  0 is  the  angle  with  the  x.  axis  of  the  projection  of  p on  the 
Li  L >v 

(x  ,x  )-plane.  It  is  noted  that  the  definition  of  p(s)  given  by  Eq.(l7) 

L Z.  *v 

is  completely  analogous  to  Eq.(3),  with  the  understanding  that  all 
quantities  now  are  functions  of  the  arc  length  s.  Considering  some 
point  s on  the  edge,  it  follows  that  the  phase  at  point  s can  be 
expressed  as 

, s 

S (s)  s S(0,0,s)  = S(0,0,s)  + cos[0  (s)lds 
O Li 

S 

For  small  x^  and  x^ , the  geometrical  optics  part  of  the  diffracted 
field  corresponding  to  the  inciil^fit  field  given  by  Eq.  (16)  can  now 
readily  be  constructed.  We  find 

~go(xl’x2’s)  = Ao(s)  exp[ikLSo(s)l  Wg0(xi>x2)  + °(1) 
where  is  defined  by  Eq.(5).  The  o(l)  term  in  Eq.(19)  contains 

the  effect  of  curvature  of  the  wavefront  at  the  edge.  Also 
A^Cs)  = A(0,0,s) 

The  expression  given  by  Eq. (19)  is  taken  as  the  outer  solution  in  a 
matching  of  two  solutions. 

Examination  of  the  geometrical  optics  part  of  the  displacement  field, 
Eq.(19),  motivates  the  following  choice  of  the  displacement  in  a small 


(16) 


(17) 


(18) 


(19) 


(20) 


region  near  the  edge  (the  inner  solution) 
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~is(yl’y2’S)  = exP[ikLsQ(s)  V(yl,y2’s)  + °(kL1)^ 
where  the  inner  variables  >\  are  given  by 

yi = Vi  > 1 = 1>2 

The  choice  of  the  stretching  factor  k is  based  on  the  requirement  that 

Li 

the  s-derivatives  of  the  inner  solution  should  be  of  the  same  order  of 
magnitude  as  the  y^-der ivatives . 

In  the  coordinate  system  (y  ,y  ,s),  the  metrics  are  h = h = k 3 , 

L Z Lc.Lt 

and  h_  = 1 + 0(k  X).  We  remark  that  the  curvature  of  the  edge  is  contai 

■J  Li 

in  the  0(k  3)  term  in  the  metric  h_.  With  j,k  = 1,2  and  implying  the 
L J 

summation  convention . the  elastodynamic  displacement  equations  of  motion  to 

8 o 
2 r 


(21) 


(22) 


ned 


leading  order  in  k become  upon  substitution  of  (21) 

Li 


kI  V,kk+(cL/cT)2(vk,ki+  1 V3,j  C0S\  + V ‘ Vj  CPS\^(kL)  = ° (23) 


kL  V,kk+(cL/cT)2(i  \,k  COSV  V3  si°2V  - V3  cos\]+0(kL)  = 0 (24) 

On  the  surfaces  of  the  crack,  the  conditions  of  vanishing  surface  tractions 
yield 

kL  [X  (\,k  + 1 V3  COSV  6j2  + U (vj , 2 + V2,j^  + °(1)  = ° 

kL  (v3  2 + L v2  COS0L^]  + CK1)  = 0 

Equation  (21)  defines  the  inner  solution 


(25) 

(26) 


The  matching  of  the  inner  and  outer  solutions  is  carried  out  on  an 

-3/4 

overlap  region  for  which  x.  = k e.  and  y.  = lOg.,  where  £.  satisfy 

1 Li  1 L Li  L 1 

2 2 k 

the  requirement  0 < (e^+e,)"1  < *>  . The  coefficients  -3/4  is  fairly 
arbitrary,  except  that  its  absolute  value  should  be  between  0.5  and  1.  The 


matching  condition  is 

-3/4  , -3/4  , ..  \ , \ . 

J*go(kL  *1*  kL  s2  ’ S _ ~is(kL®l,kLe2,S 


.(D 


(27) 


where  u and  u.  are  defined  by  Eqs.(19)  and  (21),  respectively. 

~go  ~LS 


Equation  (27)  holds  in  the  overlap  regions,  except  on  the  boundary  of  the 
shadow  region,  and  the  boundary  of  the  region  of  reflected  rays.  We 
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remark  that  the  diffracted  field,  which  is  discussed  in  detail  in  Ref. 16  , 
enters  the  matching  condition  in  subsequent  terms,  but  not  in  the  terms 
of  the  order  that  are  considered  here. 

Let  us  now  return  to  the  diffraction  of  a plane  wave  of  the  form 
given  by  Eq.(2),  by  a semi-infinite  crack.  The  total  displacement  field 
for  this  problem  depends  on  only  through  the  factor  exp (ik^p^x^) . Thus 
this  displacement  field  may  be  represented  by  an  expression  of  the  form 
u(x1,x2,x3)  = exp(ikLp3x3)  W^.x^ 

It  can  now  be  verified  that 

v(yi,y2,s)  = Aq(s)  W(yi/kL,y2/kL) 

satisfies  the  displacement  equations  of  motion  (23)  and  (24)  as  well  as 
the  boundary  conditions  (25)  and  (26)  exactly  to  the  order  indicated. 
Moreover,  with  Aq(s)  defined  by  Eq.(20),  it  is  found  that  the  matching 
condition  (27)  is  satisfied  to  the  order  indicated,  since  for  y.  = k^e-j 

1 L 1 

W(k  ^//^e.,k  ^^6  ) asymptotically  equals  W (kT^^en,  k'^^e,,).  On  the 
/v  L 1 L Z L 1 L Z 

basis  of  these  observations,  Eq. (21)  then  yields  for  the  inner  solution 

U.  (x,,  x9ls)  ~ A (s ) expfik  S (s)]  W(x  ,x  ) + 0(k'1) 

“■is  i l o L o ~ 1 z L 

This  result  implies  that  the  stress  intensity  factor  at  a point  s on  the 

edge  is  equal  to  the  field  of  the  incident  ray  at  this  point,  times  the 

stress  intensity  factor  divided  by  the  phase  exp(ik^p3x3)  for  a plane  wave 


(28) 


(29) 


(30) 


of  unit  amplitude  which  is  incident  with  the  same  propagation  vector  p. 

3.  Elastodynamic  Stress  Intensity  Factors 
Elastodynamic  stress  intensity  factors  for  a straight-edged  semi- 
infinite crack  were  analyzed  by  Achenbach  and  Gautesen  [7^.  The  diffraction 
of  a plane  longitudinal  wave  of  the  form 
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u X>  = p exp(ik  p.x)  ( 

which  is  incident  under  an  arbitrary  angle  with  the  edge  of  the  crack, 
generates  Mode  I,  II  and  III  stress  intensity  factors.  An  incident  trans- 
verse wave  of  the  form 

u = d exp(ikTp.x)  ( 

generally  also  generates  stress  intensity  factors  in  all  three  fracture 
modes.  In  Eq.  (32)  the  index  m can  assume  the  values  m = 2 and  m = 3, 
which  refer,  respectively,  to  vertical  and  horizontal  polarization  of  the 
incident  wave,  relative  to  the  plane  spanned  by  i^  and  p,  see  Fig.l.  Thus, 
d<2)  = d<3)  vp  ~ C 

~0)  = E * ^2  7 1 * i2  I C 

For  each  incident  field  u^"*  and  u7m\  m = 2,3,  the  three  stress  intensity 

/ \ 

factors  may  be  represented  by  a vector  KF  . For  an  incident  longitudinal 
wave  the  components  of  K^^are  presented  in  Ref.  [7].  The  stress  intensity 
factors  generated  by  incident  transverse  waves  can  be  computed  by  the 
methods  of  Ref. (7].  For  convenience,  the  components  of 

K(l)  _ K(!)  exp (-ik  p.x.)  (! 

^ ^ O L J J 


„(m)  v (m)  . , 

K = K exp (- ik  p.x.) 

~ ~ O I J J 


m = 2,3 


are  listed  in  Appendix  A.  Here  are  the  Mode  I,  II  and 

'•l  ~L  •'-'J 

III  stress  intensity  factors,  respectively. 

In  this  paper  we  are  interested  in  the  stress  intensity  factors 
generated  at  the  edge  of  a curved  crack,  see  Fig. 2 by  an  incident  wave 
of  the  form 

u = A(x1,x2,s)  e^m^  (s)  exp[ik('m')S(x1,x2,s)] 
where  m = 1 signifies  a longitudinal  wave,  when 
e(1)  P = VS  , k<X>  = kT 


while  m = 2 signifies  a transverse  wave  with  vertical  polarization 


r 


,(-) 


= k. 


O') 


• uul  m 3 refers  t » a transverse  wave  with  horizont  ! n . 1 a i : ■■ 


U’  - ,0>  l.(3)  . 

«'  = d , k = k 


r 

. f v.  shown  in  the  previous  section  that  for  an  incident  wave  oi 
arbitrary  shape  the  stress  intensity  factor  at  a p ont  on  crack 
edge  is  equal  t the  field  of  the  incident  ray  at  this  point  multi- 
plied by  the  stress  intensity  factor  for  a corresponding  plane  wave 
divided  by  the  phase  exp(ik^m^p^x^) . Thus 

K = ft  A e exp(ik^  'S  ) 

~ 0 •“O  o 

where  A e*'™'1  expfik^'s  ) is  the  field  of  the  incident  ray  at  a point 
o ■»  o o 

on  the  edge,  and  ),  is  a matrix,  whose  components  are 

K.k-^\uW2)  +K<3V3> 

jk  j rk  j k J k 

The  coefficients  are  defined  by  Eqs.(35)  and  (36). 

Several  specific  examples  of  elastodynamic  stress  intensity  factors 
generated  by  body  waves  are  discussed  in  Section  5. 

4.  Surface  Waves 

The  results  of  the  preceding  sections  can  easily  be  extended  to  the 
computation  of  elastodynamic  stress  intensity  factors  which  are  generated 
at  the  diffraction  of  Rayleigh  surface  waves  by  a crack  edge.  In  this 
paper  we  consider  stress  intensity  factors  due  to  surface  waves  which  are 
incident  on  an  edge  via  both  faces  of  a crack.  Then,  it  is  convenient 
to  distinguish  between  symmetric  and  antisymmetric  motions  relative  to 
the  plane  of  the  crack.  Only  the  case  of  symmetric  motions  is  considered 
in  detail,  but  the  corresponding  results  for  antisymmetric  motions  can 
be  obtained  analogously. 

Following  Eqs.(3.12)  and  (3.13)  of  Ref.  [8]  with  appropriate  minor 
adjustments  in  notation,  the  time-reduced  displacements  for  the  incident 


f 4< 


(41) 


(42) 
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surface  wave  may  be  written  in  the  forms 

R _ r / 1 9 2/  -\  ~*a|x  | , 2.  2.  -U)b|x  , . 

u = f (l-2s  / s ) e 2 + 2(s  /s)  e 2 ] V0sgn(xo) 


R 1 r,0  2 2 -uoa  x . 2.  -li'b  x 

u.  = “[l2^-^)  e 2 - 2 (ab/ s ) e 2 

X cl 


R T 


R T 

where  j = 1,3,  and 

V = expfiaT-s  (x  sinG.+  x cosg. ) ]} 

Z K 1 1 j 1 


2 

h 


, 2 2,% 

SR_SL 


2 2 % 
b = (sR-sT)^ 


Thus,  in  the  geometry  shown  in  Fig.l,  the  surface  wave  is  incident 
on  the  edge  under  an  angle  0^with  the  edge. 

The  diffraction  of  surface  waves  has  been  discussed  in  some  detail 
in  Refs.^8J  and  [9],  If  we  denote  the  time  reduced  displacement  in  the 
x^-direction  on  the  faces  of  the  crack  by  U^Cx^x,,),  the  exponential 
Fourier  transform  with  respect  to  x^  of  is  given  by  Eq.  (44)  of 


Ref.(8j.  Here  we  slightly  rewrite  this  expression  in  the  form 

, . , , E(s„sing.) 

1 1 1 1 l R l -lUI-'x.. 


*"2^^  iU)  .F-s  sin0.+  iO  F+  s sri0.+  iO, 


R 


R 


E(f) 


In  Eq. (47) 


1 = SRCOS0i 

E(|)  = K+(F,r)/r(sR-T2)^+5^ 

where  k”*"  ( ^ , T) ) is  defined  by  Eq.(4.6)  of  Ref.  [8J.  The  incident  and 

reflected  surface  waves  correspond  to  the  contributions  of  the  poles 

in  inversion  integrals.  Thus,  the  contribution  at  F = s sin0.-iO  just 

reproduces  the  incident  wave  given  by  Eq.(43),  while  the  pole  at 

F = -s  sin0.-iO  gives  the  reflected  surface  wave.  It  is  easily  veri- 
K 1 

fied  that  tie  reflection  coefficient  can  be  written  as 


(43) 

(44) 

(45) 
(46a , b ) 


(47) 

(48) 

(49) 


^ * -E(sRsin0.  J/EC-s^inep 

By  virtue  of  asymptotically  valid  relations  between  a Fourier 
transform  as  | gt' | -»  <E  and  the  inverse  transform  as  x]  — 0,  the 
following  expression  for  the  Mode  I stress  intensity  factor  has  been 


derived  as  Eq.(66)  of  Ref. [7]: 


„R  ,3/2  3/2  2,  2,  i.rr/4,  . -+,  . ... 

K0  = 2 u)  u(l-s  /s  ) e lim  U2(C)  (5. 

| - 00 

R R 

By  substituting  Eq.(47)  into  Eq.(51),  we  obtain  = K exp (- iX'x^ ) , where 
„R  „ 5 / 2 \ 2,  2.  . . . . -in/4 

K = -2  cju  h.  (l-sL/sT)  s^inG^  E(s  smG^  e (5. 

An  extension  to  incident  surface  waves  with  curved  wavefronts, 
which  are  incident  on  curved  crack  edges,  is  now  completely  analogous 
to  the  case  of  body  waves.  Suppose  the  time  reduced  incident  field  is 
of  the  form 

u2  = ±A(x1,s)  exp[-i!«sR  SCx^s)]  (5: 

where  the  displacement  is  positive  on  the  upper  face.  The  result  then 
simply  states  that  the  stress  intensity  factor  at  a point  on  the  edge 


KT  = A K expl-ia’  s„S  1 
I o L R o 


where 

A = A(0,s)  , S = S (0, s) 

o o 

Surface  waves  of  the  type  given  by  Eq.  (53)  can  be  generated  by  direct 
application  of  disturbances  to  the  faces  of  the  crack,  as  well  as  by 
diffraction  processes.  Examples  of  both  cases  are  given  in  the  next 
section. 


5 . Examples 


In  this  section  several  applications  of  Eqs.(41)  and  (54)  are  pre- 
sented. We  will  consider  examples  with  curved  as  well  as  with  plane 
wavefronts,  and  with  curved  as  well  as  with  straight  crack  edges.  Both 
steady-state  time  harmonic  and  impact  excitations  will  be  considered. 
Anti-Plane  Motions 

The  first  example  is  concerned  with  the  stress  intensity  factor  at  the 
tip  of  a semi-infinite  crack  due  to  an  anti-plane  line  load  of  P force 
units  per  unit  length  and  time  variation  exp(-iU)t),  which  is  applied 
parallel  to  the  crack  as  shown  in  Fig. 3.  This  elastodynamic  problem  is 
two-dimensional,  and  the  time  reduced  anti-plane  displacement  u^Cx^.x^) 
is  governed  by 

V2u3  + k2u3  = - (P/u)  6(r  - rp)  (56) 


The  solution  to  Eq.(56)  is 

"3(I-V  -U  “o'Xlt-  Jpl> 

For  large  frequencies  we  may  write 


u3(r>rp) 


i P 2 


l h 


4 * -nkT  lr£P!' 


ITT 

A 


(57) 


(58) 


exp[ikTlr-rp| 

Equation  (58)  represents  the  incident  field. 

For  an  incident  plane  horizontally  polarized  transverse  wave  of  unit 
amplitude,  the  elastodynamic  stress  intensity  factor  was  presented  by 
Eq.(90)  of  Ref. [7]  as 

t _ I'  T-T  / /, 

(59) 


(*?>),  - *** 


where  g^,  is  the  angle  of  incidence.  For  the  problem  at  hand  we  have 

p = cosg  = r(x  ) - (x.)n'/ !r  -r  I.  Equation  (59)  can  then  be  re- 

1 1 i U i F *o  ~P 


written  as 


(60) 


-O)' 


= p.  (2k  ) (1  - Pl)2  e 


\ -in/4 


v o 3 " '*  rV 

Equation  (41),  together  with  Eqs.(58)  and  (60)  then  yields  for  the  line 

load  problem 

K3  = 2 P T7'-r”|)%(1"Pl)^eXp[ikT  IEo'^P  1 ] (61) 

'*o  ~p  ■ 

This  result  agrees  with  the  near-tip  field  computed  from  the  asymptotic 
expansion  of  the  exact  solution  of  an  analogous  problem,  see  Ref. f 101, 
p.  329,  Eq. (8.74). 


For  the  second  example  we  consider  a crack  of  finite  length  U , 
subjected  to  a plane  horizontally  polarized  transverse  wave  of  the  form 

(u3)inc  = A exp[ia'(s  x cos0  + s^sine^,  - t)  ] (62) 

The  geometry  is  shown  in  Fig. 4.  To  compute  the  Mode  III  elasto- 
dynamic  stress  intensity  factor  at  high  values  of  we  trace  an 

incident  ray  as  it  undergoes  diffractions  at  the  two  crack  tips.  The 
primary  diffraction  produces  diffracted  rays  on  the  faces  of  the  crack, 
whose  fields  for  large  frequency  and/or  distance  are  easily  computed,  as 
shown  in  Ref. [6".  These  fields  are  subsequently  diffracted  at  the  other 
crack  tip,  thus  producing  another  set  of  rays  on  the  crack  faces.  This 
process  continues,  and  one  can  visualize  systems  of  rays  travelling  back 
and  forth  between  the  two  crack  tips.  Since  the  fields  decay  as  (ous^p)  2, 
it  is  only  necessary  to  consider  a few  secondary  diffractions.  In  fact, 
in  the  computations  presented  here,  only  one  secondary  diffraction  is 
inc luded . 

The  primary  diffracteu  field  of  the  incident  ray  at  0^,  see  Fig. 4, 
is  the  diffracted  field  for  a semi-infinite  crack.  This  field  was  dis- 
cussed in  Ref. [6],  On  the  upper  surface  of  the  crack  (i.e.,  at  0 = 0) 
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the  exponential  Fourier  transform  of  the  anti-plane  displacement  is 
given  by  Eq. (89)  of  Ref. [7]  as 


U3  = (2krf)^  A ' iX3^2(P+srrcoseT)  (|+sT)'2  1sin(0T/2) 


(63) 


Analogously  to  the  results  given  by  Eqs . (65) , (70)  and  (71)  of  Ref. r6’, 
the  far  field  (or  high  frequency)  inversion  at  0 = 0 can  be  computed 


03- 


where  the  diffraction  coefficient  D^,  is 

n+  - __L_  % 1 i"/4 

T ~ v2rwsT  sin (0/2) 


0 < 0^  < 


(64) 


(65) 


The  primary  diffracted  field  at  0?  is  obtained  by  taking  into  account 
that  the  angle  of  incidence  measured  from  the  shadow  side  of  the  crack 


is 


as  shown  in  Fig. 4.  For  the  primary  diffraction  we  can  then 


write 


U3  = ("x]+2,l->  “ DT  A e 


iaJ(-sTx1+2sT£-t) 


where 


1 '% 


in/4 


T \2t«jst  cos(0T/2) 


0 < eT  < tt 


(66) 


(67) 


The  corresponding  displacements  on  the  lower  crack  faces  are  equal  in 
magnitude  but  opposite  in  sign  to  Eqs. (64)  and  (66). 

For  the  primary  incident  ray  the  stress  intensity  factor  at  = 0 
is  given  by  Eq.(59),  which  is  valid  in  the  range  0 < <.  n.  For 

0^,  = 77  , the  waves  approach  the  crack  tip  from  the  side  of  the 
crack  face.  Then,  not  only  the  incident  wave  but  also  the  reflected  wave 
is  incident  on  the  crack  tip.  For  this  system  of  incident  and  reflected 
waves,  the  displacements  on  the  crack  face  are  actually  twice  as  large 


as  would  correspond  to  the  incident  wave  only.  This  observation  implies 
that  for  an  incident  wave  which  is  given  as  a surface  displacement  on 
the  lower  face  of  the  crack,  the  elastodynamic  stress  intensity  factor 


is  half  the  value  given  by  Eq.  (59)  for  0^,  = tt. 

On  the  basis  of  the  foregoing  observation,  the  stress  intensity 
factor  at  x^  = 0 due  to  the  primary  diffracted  wave  generated  on  the 
lower  crack  face  (which  is  minus  the  expression  given  by  Eq.  (66) ) , then 
follows  by  setting  © = tt  in  Eq. (59)  and  substituting  the  appropriate 

amplitude,  as 


Km  = - Hk^  e_1"/4(2£)'^  D~  A e2lU,S/ 

Multiplying  this  result  by  2,  to  account  for  the  diffracted  wave  on 
the  upper  face,  which  is  given  by  Eq.(66),  and  adding  the  result  to  the 
primary  stress  intensity  factor  given  by  Eq.(59),  we  obtain  at 


X!  = ° 

|Kin|~  2uAk^!e'iTT/4[sin(0T/2)  - (21  )"%d‘  e2iU)s/  ] | 

It  is  easily  shown  that  the  incident  displacement  wave  given  by 
Eq.  (62)  corresponds  to  a stress  wave  of  magnitude  t , where  = uAk^. 

Let  us  now  consider  the  special  case  of  normal  incidence,  when  0^  = tt/2 
Then  we  have 


' 2 -iTT/4  L.\%  i(2k  £+rr/4)_ 

vky  ,e  -'ny/  e T 


which  can  be  further  simplified  to 

Wv%l  ~ - y ^os(2kT«  +W4)] 

This  stress  intensity  factor  was  computed  by  other  methods  by  Mai  [11]. 
The  results  are  compared  in  Fig. 5.  It  is  noted  that  the  agreement  is 


(68) 


(68) 


(69) 


(70) 
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excellent  for  k^f,  > 1.5. 

In-Plane  Motions 

As  an  example  of  the  much  more  complicated  class  of  in-plane 
problems  we  will  consider  the  elas todynamic  stress  intensity  factor 
generated  by  the  diffraction  of  a longitudinal  wave  by  a crack  of 
length  2£. 

Similar  to  the  anti-plane  problem  the  e lastodynamic  stress  intensity 
factor  is  computed  by  tracing  rays  as  they  are  diffracted  by  the  two 
crack  tips.  There  are,  however,  some  interesting  and  important  differences 
with  the  case  of  anti-plane  motions.  The  most  important  difference  is 
related  to  the  significance  of  Rayleigh  surface  waves  in  the  primary  and 
secondary  diffractions.  The  primary  diffraction  generates  not  only  rays 
of  longitudinal  and  transverse  motion,  but  also  surface  motions  near  the 
faces  of  the  crack.  It  is  a general  feature  of  in-plane  problems  that 

-h 

bulk  waves  decay  as  (distance)  , while  surface  waves  do  not  decay  at 
all.  Thus,  on  the  faces  of  the  crack,  the  surface  motions  tend  to  pre- 
dominate the  longitudinal  and  transverse  motions  due  to  the  body  waves. 
This  means  that  it  is  necessary  to  consider  stress  intensity  factors  due 
to  surface  waves,  which  were  discussed  in  Section  4.  Since  there  is  no 
spatial  decay,  except  for  the  fact  that  the  reflection  coefficient 
given  by  Eq.(50)  is  smaller  than  unity,  it  is  necessary  to  sum  the 
influence  of  all  secondary  diffractions  of  the  surfaces  waves.  Here 
we  present  the  details  for  the  case  of  normal  incidence. 

The  Rayleigh  surface  wave  generated  by  the  primary  diffraction  of 
the  incident  ray  follows  from  Eq.(46)  of  Ref. [6].  For  the  problem  at 


hand  we  have 
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(u2>inc  = A exPriu(sLx2-t)] 

Thus  = P3  0 , p9  = 1.  The  exponential  Fourier  transform  of  the 
surface  displacement  for  x^  a 0 then  follows  from  Eq.(46)  of  Ref.^"1 
as 


iST/sL 


2^J(s^-s^)%  5 


(SL+^' 
SR+  5 


A 

K+(?) 


(71) 


(72) 


The  surface  wave  contribution  comes  from  the  pole  at  5 = - s 


R 


We  find 


R _ . i(«(s  x - tut) 

u2  = DR  A e R 1 


(73) 


where 


n ■ !l  SR/SL-1  1*  1 

°R  1 s_  „ „ 2 , 2S  ..+ 


(74) 


R L2(l-sL/sT)  K (-sR) 


The  Mode  I stress  intensity  factor  due  to  the  primary  diffraction 
immediately  follows  from  Eq.(41)  by  setting  p^  = p^  = 0 , p^  = 1.  We 
find 

K®  = K^1)  A (75) 

where  can  be  obtained  either  as  a special  case  of  Eq.(A.2),  or 

by  employing  Eq.(51)  together  with  Eq.(72).  The  result  is 

K^1}  = 2uA(s^/sJ-l)%(uusL)%  e"W4  (76) 

The  secondary  diffractions  at  point  0^  are  produced  by  Rayleigh  sur- 
face waves.  The  phases  of  subsequent  contributions  as  a ray  travels 
back  and  forth  between  the  crack  tips  increase  with  2U)s  f with  each 
diffraction,  i.e.,  they  are  exp(2ious  £),  exp(4iu>s  ),  etc.  The 
magnitudes  of  subsequent  contributions  are  reduced  by  multiplication 
by  the  reflection  factor  R^  given  by  Eq.(50).  Thus,  for  the  contribu- 
tions of  the  Rayleigh  waves  to  the  stress  intensity  factors,  we  have 


KR  = K^A,  where 


KR  = Dr  KRre2iU)S/+RRe4iU,SR£-H^e6iU,V+...: 
where  is  defined  by  Eq.(52).  Equation  can  be  written  in  closed 


(77) 


form  as 


R , Dr  Kq  exp(?.iu»sR0 
L-R  exp(2iats  ■ ) 

K K 

The  total  stress  intensity  factor  is 

*!-«?♦  >4* 


(78) 


(79) 


2 2 

Numerical  results  were  worked  out  for  s“ / s ^ = 3,  which  corresponds 

to  a material  with  a Poisson's  retain  v = 0.25.  For  that  case  we  have 

K+(s  ) = 0.91116  and  K+(-s  ) = 1.9009.  The  absolute  value  of  KT  is 
K R 1 

compared  in  Fig. 5 with  numerical  results  obtained  by  Hal  [111.  It 
is  noted  that  the  agreement  is  very  good  for  > 1.5. 

Normal  Point  Loa 4s  on  a Semi-Infinite  Crack 

As  an  example  of  the  computation  of  elastodynamic  stress  intensity 
factors  generated  by  surface  waves  only,  we  consider  the  case  that  the 
faces  of  a semi-infinite  crack  are  subjected  to  equal  and  opposite  time- 


harmonic  point  loads  applied  at  = ± 0 , and  = 0 , see 

Fig. 6.  The  boundary  conditions  at  x^  = 0 then  are 
T 22  = ~ N 6(Xj-xo)  5(x.j)  exp(-iait) 

0 


21 


23 


(80) 

(81) 


Motions  generated  by  a normal  point  load  on  a half-space  have  been  studied 

2 2 2 

in  detail.  It  is  well  known  that  for  »s  r > > 1,  where  r = (x.-x  ) +x_, 

R I o J 

the  predominant  component  of  the  surface  displacement  is  the  Rayleigh 
surface  wave.  The  normal  surface  displacement  is  given  by  Eq. (2-120)  of 
Ref. [12].  Taking  into  account  that  in  the  present  paper  we  use  a 


negative  time  exponent,  exp(-iu)t),  the  time-reduced  displacement  can  be 
written  as 


u2  „ A<r)  e‘<"V+  "/4)  (82) 

where 

s^/s2  2.% 

A(r)  R L. V.  ,R  y (83) 

A(  ’ U \2vr,'  R '(s  ) (83) 

R 

In  Eq.(83),  R(c)  is  the  Rayleigh  function 

R(r)  = <2C2-s^)2  +4C2(s^-  C2)^(s2-  C2)^  (84) 

and  the  prime  denotes  the  derivative  with  respect  to  C. 

In  the  context  of  the  discussion  of  Section  4,  and  referring  to 

Fig.  6,  we  can  now  define  at  a point  x^  on  the  edge 

2 2 \ 

sin0  = x /r  ; cosp  = x„/r  ; r = (x  + x_K  (85) 

10  10  O J 

By  comparing  Eqs.  (53)  and  (82)  we  can  now  immediately  identify  A(x^,s) 
and  S(x  ,s).  By  virtue  of  Eq.  (54)  the  stress  intensity  factor  may 
then  be  written 

K = A(r)  KR  exp[ix(s  r + tt/4)  ] (86) 


This  stress  intensity  factor  was  also  computed  in  Ref. [71  by  a direct 
asymptotic  evaluation  of  integrals  in  an  expression  for  the  exact 
solution.  It  can  be  shown  that  Eq,  (86)  agrees  with  Eq. (104)  of  Ref.  [7], 
Pennv-Shaped  Cracks 

For  the  final  example  we  consider  a penny-shaped  crack  under  the 
influence  of  an  axially  symmetric  torsional  wave.  In  cylindrical  co- 
ordinates (r.gjX^),  axially  symmetric  torsional  wave  motions  are  defined 

by  a single  displacement  component  u (r,9,x„)  which  is  independent  of  the 

0 ^ 

angular  variable  p.  The  only  non-zero  stress  components  are  t and  t 

20  rfl 

In  the  problem  considered  here  the  axis  of  symmetry  is  normal  to  the  plane 


of  the  crack,  and  passes  through  the  center  of  the  penny-shaped  crack. 


The  geometry  is  shown  in  Fig. 7.  The  incident  stress  wave  is  of  the  form 
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V> 


Xi 


?\X- 


— H(cTt  - x2) 


(87) 


Thus,  this  is  a transient  stress  wave,  which  strikes  the  crack  at  time 
t = 0.  Writing  the  near-tip  stress  field  in  the  form 

, . 1 


K3(t) 


(88) 


2e  (2n)^  (r-a)^ 

the  first  term  in  an  expansion  over  time  follows  from  the  results  of 

Kennedy  and  Achenbach  F13,  £q.(3.11)]  as 

K3(t)  = 2To(2cTt/n)^  [l  + 0(t)]  (89) 

Now  we  will  approach  this  problem  with  the  theory  presented  in  this 

paper.  The  incident  wave  (u  ) . can  be  expressed  as 

0 me  r 


(u0)inc 


CTTor 


j — exPriu)(x2/cT-t)]duJ 


(90) 


~ no  (JJ 


Treating  the  quantity  (c^t^r /2TT  ViatO  )exp(iXx2/c^,)  as  the  incident  wave, 
it  is  observed  that  this  wave  is  equivalent  to  an  incident  plane  hori- 
zontally polarized  time-harmonic  wave,  which  is  normally  incident  on  the 
crack,  i.e.,  p^  = 0.  It  then  follows  by  employing  T5qs.(41)  and  (60)  that 


/2u>'*  -rT/4  CTTo 

K„  = U,  — e r 

3 'V  2n  vm2 


(91) 


Now  we  can  introduce  the  time  factor  exp(-iU)t)  in  Eq.  (91)  and  substitute 


the  result  in  the  Fourier  superposition  integral.  Since 


- 3/2  -ixt 


dt  = 4(tTt)  e 


h 1^/4 


(92) 


the  evaluation  of  the  superposition  integral  yields  Eq.(89),  and  the 
results  obtained  by  the  method  of  this  paper  agree  with  those  of  Ref. [131. 
This  example  also  confirms  the  expected  result  that  the  large  frequency 


\ 

i 
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approximation  in  the  frequency  domain  leads  to  a small  time  approximation 
in  the  time  domain.  In  this  context  small  time  means  small  after  arrival 
of  the  wavefront  which  separates  the  disturbed  and  the  undisturbed  regions 
of  the  solid. 
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APPENDIX  A 

ELASTODYNAMIC  STRESS  INTENSITY  FACTORS  FOR  THE  INCIDENCE 
OF  PLANE  WAVES  ON  A SEMI- INFINITE  CRACK 

In  this  Appendix  we  have  listed  the  components  of  ^ and 

m = 2,3,  which  are  defined  by  Eqs.  (35)  and  (36).  For  economy  of  pre- 
sentation the  slownesses  s , n = 1,2,3,  are  introduced  as 

s L = l/cL  , s9  = l/cT  and  s3  = 1/c  , 

where  cR  is  the  velocity  of  Rayleigh  waves.  In  the  following  ex- 
pressions the  index  (j)  can  assume  the  values  (1)  (incident  longitudinal 
wave)  and  (2)  (incident  vertically  polarized  transverse  wave).  The 
index  (3)  refers  to  an  incident  horizontally  polarized  transverse  wave. 
The  stress  intensity  factors  are: 
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